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ABSTRACT 


An  investigation  was  conducted  to  find  suitable  devices  and  techniques 
for  measuring  heat  transfer  rates  in  the  base  areas  of  rocket -powered 
missiles.  A  slug  mass  calorimeter  was  selected  to  measure  total  (convec¬ 
tive  and  radiant)  heat  rates,  and  a  black-body  cavity  was  selected  to  meas¬ 
ure  radiant  heat  rates.  Similar  instruments  have  been  used  in  past  test 
programs;  however,  the  simplifying  methods  of  analysis  as  used  in  the  past 
with  these  instruments  may  result  in  considerable  errors.  The  approach 
taken  in  this  report  was  to  determine  the  source  heat  rate  by  writing  a  gene¬ 
ral  heat  balance  equation  as  a  function  of  time  about  the  heat  sensing  unit. 

An  experimental  program,  utilizing  a  box-type  muffle  furnace  to  achieve  a 
known  source  heat  rate,  was  conducted  to  determine  the  accuracy  of  the 
mathematical  analysis  used  herein.  The  results  of  the  experimental  pro¬ 
gram  illustrated  that  heat  transfer  rates  may  be  measured  with  reasonable 
accuracy  providing  basic  fundamentals  are  adhered  to  and  good  engineering 
practices  are  utilized  in  the  instrument  design.  In  addition,  the  instruments 
and  techniques  as  presented  in  this  report  should,  in  general,  apply  to  other 
testing  applications. 
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NOMENCLATURE 

A  Heat  transfer  area,  ft2 

cp  Specific  heat,  Btu/lbm-deg 

D  Diameter,  ft 

F  Angle  factor 

h  Film  coefficient,  Btu/ft2-hr-deg 

k  Thermal  conductivity,  Btu/ft2-hr-deg/ft 

JL  Height  of  black-body  cavity  cone,  ft 

m  Mass,  lbm 

Q  Heat  transfer  rate  at  any  time,  units  as  specified 

r  Radius,  ft 

S  Internal  surface  area  of  the  black-body  cavity,  ft2 

S0  Equivalent  surface  area  of  a  sphere  of  the  same  depth 

as  the  black- body  in  the  normal  direction  of  the 
aperture,  ft2 

s  Cross-sectional  area  of  the  aperture,  ft2 

T  Temperature,  °R 

T'  Temperature  after  a  small  time  increment,  °R 

t  Thickness,  ft 

y  Distance,  ft 

a  Thermal  diffusivity  =  -p  cp"  *  ft2 /hr 

/3  Time  as  used  in  Eq.  (3),  hr 

S  Increment  of  distance,  ft 

«  Emissivity 

Apparent  emissivity  of  a  black-body 
9  Time,  hr 

A  Integration  variable  (time) 

p  Density,  lbm/ft3 

a  Stephan- Boltzmann  constant,  0.  1713  x  10"8 

Btu/ft2-hr-°R4 
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<f>  Surface  temperature  function 

x{t  Half  angle  of  black-body  cone 

SUBSCRIPTS 

l,  2  ... .  Refer  to  location  or  component 
A  Calorimeter 

b  Insulator 

c  Thermocouple 

D  Heat  sensor 

g  Gas 

i  Initial 
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1.0  INTRODUCTION 


Attempts  to  measure  heating  rates  during  various  simulated  alti¬ 
tude  tests  of  missile  bases  at  AEDC  and  elsewhere  have  not  been 
entirely  successful,  primarily  because  of  inadequacies  in  heat  trans¬ 
fer  measuring  equipment  and  techniques.  The  specific  purpose  of 
this  program  was  to  produce  suitable  devices  and  techniques  for  meas¬ 
uring  heat  transfer  rates  in  the  base  area  of  rocket-powered  missiles; 
however,  the  content  of  this  report  is  not  restricted  to  this  application 
and  may  find  use  in  other  types  of  test  programs,  such  as  nose  cones, 
re-entry  bodies,  and  duct  flow.  This  investigation  was  conducted  in 
the  Rocket  Test  Facility  (RTF),  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC). 

The  intent  of  a  given  test  program  will  determine  the  specific 
parameters  which  are  to  be  measured,  and  for  the  purpose  of  this  re¬ 
port,  it  was  assumed  the  following  missile  base  parameters  are  desired: 

(1)  convective  heat  rates,  (2)  radiant  heat  rates,  and  (3)  base  gas  film 
coefficients.  Film  coefficients  are  particularly  important  if  experi¬ 
mental  data  are  to  be  presented  in  a  general  form  to  be  useful  in  the 
formulation  of  theoretical  analysis.  Devices  and  methods  of  determining 
the  above  parameters  are  varied  and  numerous,  but  conditions  which 
are  associated  with  the  measurement  of  heat  rates  on  a  missile  base 
limits  the  selection  of  appropriate  instrumentation.  These  limitations 
dictate  that: 

1.  The  instruments  must  be  of  a  flight- weight  type  and  robust  to 
survive  rough  handling  and  the  intense  vibrations  which  are 
associated  with  the  start  transients  of  rockets, 

2.  The  instruments  must  not  be  adversely  affected  by  deposition  of 
combustion  exhaust  products  on  their  surfaces, 

3.  The  instrument  used  to  measure  convective  heat  rates  must  not 
appreciably  disturb  the  thermal  characteristics  of  the  base,  and 

4.  The  instruments  should  be  designed  so  they  are  amenable  to 
mathematical  analysis. 

To  measure  the  total  heat  rate,  an  insulated  slug  mass  was  selected 
which  had  approximately  the  thickness  and  was  constructed  of  the  same 
material  as  the  base.  For  short  duration  tests,  it  is  necessary  to  duplicate 
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only  the  depth  of  the  material  which  will  be  affected  by  a  thermal  change 
during  the  test.  To  measure  radiative  heat  rates,  a  black-body  cavity 
similar  in  concept  to  the  cavity  utilized  by  various  organizations  as  a 
radiative  standard  (Refs.  1  and  2)  is  proposed.  Upon  determining  total 
and  radiant  heat  rates,  the  gas  film  coefficient  may  be  evaluated  by  meas 
uring  the  base  gas  temperature.  No  analysis  is  presented  concerning 
gas  temperature  measurements  since  several  adequate  probes  of  this 
nature  have  already  been  developed  and  analyzed  (see  Refs.  3,  4,  and  5). 

Similar  instrumentation  to  that  proposed  herein  has  been  used  on 
past  test  programs.  The  general  method  of  analysis,  however,  was  to 
assume  that  the  change  of  internal  energy  of  the  heat  sensing  unit  was 
equal  to  the  source  heat  rate  or  to  compare  experimental  results  with 
laboratory  calibrations  to  determine  the  source  heat  rate.  It  was  shown 
during  this  investigation  that  both  these  methods  may  be  greatly  in  error. 
The  approach  utilized  in  this  report  was  to  determine  the  source  heat 
rate  by  writing  a  general  heat  balance  equation  as  a  function  of  time 
about  the  heat  sensing  unit.  Because  it  was  impractical  to  try  to  antici¬ 
pate  every  type  configuration  and  the  appropriate  method  of  mathematical 
evaluation,  a  design  was  selected  which,  in  general,  would  apply  to  heat 
transfer  measurements  encountered  in  rocket  base  heating  tests.  How¬ 
ever,  it  should  be  stressed  that  each  system  should  be  analyzed  on  its 
own  merits  to  determine  what,  if  any,  valid  simplifying  assumptions  may 
be  made  to  assist  in  the  reduction  of  the  data. 


2.0  ANALYSIS  OF  SLUG  MASS  CALORIMETER 

In  the  earlier  phases  of  this  program,  Westkaemper  (Ref.  6)  ana¬ 
lyzed  the  factors  influencing  accuracy  when  using  a  slug  mass  calorim¬ 
eter  to  determine  total  source  heat  rates.  The  principal  conclusions  in 
Ref.  6  are: 

1.  The  effects  of  surface  temperature  mis-match  between  the  meas¬ 
uring  device  and  surrounding  base  is  a  source  for  considerable 
error.  The  magnitude  of  this  error  (as  stated  in  Ref.  6)  has 
since  been  found  to  be  somewhat  pessimistic  (Ref.  7);  however, 

a  temperature  mis-match  is  still  a  large  potential  source  of  error. 

2.  The  calorimeter  should  duplicate  the  surface  conditions  of  the 
segment  it  replaces. 

3.  The  calorimeter  mounting  insulation  may  be  a  source  of  con¬ 
siderable  conduction  losses. 
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Adherence  to  the  first  and  second  conclusions  is  particularly  important 
on  multi-jet  missile  bases,  since  at  present  no  adequate  theory  exists 
to  predict  the  magnitude  of  the  effect  caused  by  deviation  from  actual 
base  similarity  (change  of  base  thermal  characteristics).  It  is  generally 
not  possible  to  have  complete  similarity  and,  in  fact,  it  may  be  desirable 
to  accept  a  negligible  or  known  deviation  in  dimensional  similarity  which 
will  simplify  the  mathematical  analysis;  however,  because  of  the  present 
state-of-the-art  of  determining  missile  base  flow  fields,  any  deviations 
should  be  meticulously  analyzed.  For  instance,  a  common  practice  in 
rocket  measuring  techniques  is  to  calibrate  a  slug-mass  calorimeter 
with  a  radiant  heat  source  and  compare  experimental  temperature  histories 
to  the  laboratory  calibration  results  and  determine  the  experimental 
source  total  heat  rate.  For  this  method  to  be  successful,  either  the  sur¬ 
rounding  environment  of  the  calorimeter  would  have  to  be  controlled  to 
that  of  the  laboratory  calibration  or  the  calorimeter  encased  in  a  quantity 
of  insulation  which  would  isolate  the  unit  from  its  surrounding  environment. 
The  above  methods,  to  control  or  isolate  the  surrounding  environment, 
could  have  a  strong  effect  on  the  gas  film  coefficient  and,  hence,  on  the 
measured  source  heat  rate.  Also,  when  calibrating  with  a  radiant  heat 
source,  the  heat  rate  absorbed  is  a  direct  function  of  the  calorimeter 
absorptivity  which,  in  general,  is  a  very  sensitive  parameter. 


2.1  MATHEMATICAL  MODEL 

A  schematic  of  the  calorimeter  which  is  to  be  analyzed  is  shown  in 
Fig.  1.  The  following  conditions  and  assumptions  are  specified: 

1.  The  calorimeter  is  made  of  the  same  material  with  approximately 
the  same  thickness  and  has  the  same  surface  conditions  and  con¬ 
tour  as  the  segment  of  the  base  which  it  replaces. 

2.  The  thermocouple  leads  are  infinitely  long  and  perfectly  insulated. 

3.  Thermal  contact  resistances  at  interfaces  between  materials  are 
zero. 

4.  The  surface  emissivity  and  absorptivity  of  the  calorimeter  are 
known  and  equal. 

5.  The  calorimeter  and  the  base  are  respectively  Newtonian  (negli¬ 
gible  internal  thermal  resistance). 

6.  Thermal  properties  for  all  the  constituents  are  independent  of 
temperature. 

7.  The  thickness  of  the  insulation  is  small  compared  to  the  radius  of 
the  calorimeter. 
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8.  The  heat  rate  is  uniform  over  the  surface  of  the  calorimeter 
but  may  vary  with  time. 

9.  The  temperature -time  history  for  the  calorimeter  may  be  ex¬ 
pressed,  for  mathematical  convenience,  by  an  empirical 
equation  of  the  following  form: 

-  Ti  +  B  (1  -  e“z  6 )  (1) 

»  experimental  constant,  6R 

-  experimental  constant,  per  hr 


where 


Ta  (6) 


B 

Z 


All  data  from  base  heating  tests  encountered  thus  far  were  found 
to  fit  the  above  relation;  however,  any  other  suitable  relation 
may  be  used. 

For  any  given  time,  the  following  heat  balance  may  be  written  about  the 
calorimeter: 


QATet,i  -  Qai  +  Qa2  +  Qa3  +  Qa4 


(2) 


where 

QaTo1,i  =  total  heat  rate  to  the  calorimeter 
Qai  =  change  of  internal  energy  in  the  calorimeter 
Qa2  =  conduction  losses  through  the  thermocouple  leads 
Qa3  =  conduction  losses  into  the  insulator 
Qa4  =  re-radiation  losses  from  the  calorimeter 
Furthermore  (Refer  to  Fig.  2  for  the  location  of  the  x  and  y  axes): 


Total  = 

AA  hG  {tG(0)  -  [B  (1  -  e-Z«)  +  Tj]}  +  fA  Aa  0GR.dUtion 

(2a) 

Qai  = 

mA  cpA  (B  Z e  ) 

(2b) 

QA2  - 

,  a  d  T  A  ( ®)  | 

-kc*rC  . 

XC  1  xc=  0 

(2c) 

Q  A3  - 

-  kB5"  D,/“  [tb  U.y,  »>  -T  *(«],.»,  ^ 

(2d) 

/  A  y  [Tb  (x,  y,  0)  -  Ta  (0)]  x  -  dy 

QA4  = 

€A  AA  O-  [b  (1  -  e~zd)  +  Ti]4 

(2e) 
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2.2  SOLUTION  OF  THERMOCOUPLE  LOSSES 


Reference  8  presents  the  following  approximate  solution  for  one¬ 
dimensional  heat  flow  in  a  semi-infinite  slab  with  a  uniform  initial 
temperature: 


d  T 
dx 


x  =  0 


*'(0) 


+ 


3) 


— —  dxl 


(3) 


for  j8  >  0 
From  Eq.  (1): 

=  B  (1  -  e~zd) 

8)  =  B  [l  -  e-z<0-0>] 

and  letting  n  =  0  -  A 

<£U)  -  -  ft)  -  B  [l  - 


Now  substituting  the 'above  relations  into  Eq.  (3)  and  integrating  and 
substituting  the  result  into  Eq.  (2c),  the  following  solution  is  obtained: 


QA2 


kc 


V  ac/jr 


2BZj31/j  e~zd  + 


[  1  -  e-z(0-0n 


p 


B 


_  B  e-z  (0-/3) 


1  +  2Z  /S1/j 

P 


_L  Za  /S7,  +  -JL.  z*  jS*4 

3  15 


+  B 


2  Z  0 


V, 


_JL  z1  03/s  +  -5-  Z*  0^ 

3  15 


(4) 


for  £  >  0,  Z  > | v  0  <  1 


2.3  SOLUTION  OF  INSULATOR  LOSSES 

In  order  to  solve  Eq.  (2d),  a  two-dimensional  grid- network  was  con¬ 
structed  (Fig.  2)  and  solved  by  explicit  numerical  methods.  Only  the 
generalities  of  this  system  are  discussed  because  a  thorough  explanation  is 
presented  in  Ref.  9  and  also  some  discussion  is  presented  in  Ref.  6.  The 
total  number  of  grid  points  to  be  used  will  have  to  be  a  compromise  be¬ 
tween  the  type  data  reduction  program  and  the  particular  application  com¬ 
pared  to  the  degree  of  accuracy  required.  One  simple  means  to  size  the 
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network  is  to  apply  a  given  set  of  data  and  increase  the  network  and  com¬ 
pare  the  deviations  of  the  results,  or  the  analytical  method  in  Ref.  10 
may  be  used  to  estimate  the  truncation  error.  Having  the  experimental 
temperature -time  histories  for  the  boundaries  of  the  insulator,  the  grid- 
network  may  be  solved  to  determine  the  temperature  gradient  in  the 
insulator  adjacent  to  the  calorimeter.  A  typical  calculation  for  the  inner 
portion  of  the  insulator  (see  Fig.  2)  is 


(T/  -  T.)  =  M  (T.  -  T.)  +  M  (T„  -  T.)  +  -L-  (T,  -  T.)  +  -L  (T„  -  T.) 

where 


M  = 


N  = 


A. 


i 


l\6  lcB 


1  +  M1 2 3 * 


KB  rB 


2  \  da-Q 

sT~ 


A  typical  calculation  for  the  surface  of  the  insulator  is 

(T-'  -  T’)  -  +  -^T5-  <*•  -  T-> +  (T<  -  T->  +  H  Lv 1>l 


where 


QBTot«l  =  ^2  [TCm  ~  Ts  1  +  fB  ^2  QBH«di«tion  Interchange 

Since  the  solution  of  the  grid-network  includes  QBTot,i»  another  assump¬ 
tion  is  required.  If  the  gas  film  coefficient  for  the  insulator  and  the 
calorimeter  are  assumed  to  be  equal  and  if  the  base  gas  temperature 
history  is  known,  the  preceding  equations  may  be  solved  by  an  iterative 
method.  This  procedure  can  prove  to  be  quite  awkward;  however,  if  the 
calorimeter  configuration  is  changed  as  shown  in  Fig.  3a,  the  analytic 
solution  will  be  greatly  simplified.  In  analyzing  this  configuration,  the 
assumptions  and  conditions  are  the  same  as  previously  stated  with  the 
addition  that  the  lateral  heat  conduction  from  the  base  through  the  inter¬ 
connecting  strip  of  metal  to  the  calorimeter  is  negligible.  This  configu¬ 
ration  has  two  principal  advantages: 

1.  Since  the  temperature  history  is  known  for  all  the  boundaries 
of  the  insulator,  a  solution  may  be  obtained  directly  without  re¬ 
sorting  to  iteration  procedures,  and 

2.  The  interconnecting  strip  of  metal  will  duplicate  the  surface 

conditions  and  if  it  is  relatively  short  it  should  remain  in  tem¬ 

perature  equilibrium  with  the  calorimeter  and  eliminate  large 

surface  temperature  gradients. 
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2.4  VARIATION  OF  THERMAL  PROPERTIES  WITH  TEMPERATURE 


If  a  large  temperature  rise  is  encountered  by  the  calorimeter,  the 
assumption  of  constant  thermal  properties  for  the  insulator  will  prob¬ 
ably  not  be  valid.  Since  the  insulator  should  be  relatively  thin  and  since 
there  is  no  special  need  for  the  insulator  to  have  a  particularly  low 
thermal  diffusivity,  the  thermal  diffusivity  of  the  insulator  may  be  esti¬ 
mated  at  the  temperature  level  of  the  calorimeter  or  the  average 
temperature  of  the  calorimeter  and  adjacent  base,  whichever  is  more 
representative  of  the  temperature  level  in  the  insulator.  If  the  thermal 
properties -of  the  insulator  are  estimated  at  the  temperature  of  the  calo¬ 
rimeter,  the  equations  for  the  grid- network  would  be  affected  as  follows: 


«b  (Ta) 


(Ta) 

P B  (Ta)  cpb  (Ta) 


8* _ 1 

2  aB  <T'a)  1  +  M1 

2  A  daB  (TA) 

8* 

The  variation  of  the  thermal  diffusivity  as  a  function  of  the  calorimeter 
temperature  will  depend  upon  the  particular  insulating  material.  Care 
should  be  exercised  in  the  above  relations  to  ensure  that  the  time  incre¬ 
ment  does  not  become  infinitesimally  small  or  excessively  large  because 
of  computation  time  and  truncation  error,  respectively. 


A  6  - 

N  = 


2.5  CALORIMETER  TEMPERATURE  GRADIENTS 


It  may  not  always  be  possible  to  adhere  to  the  assumption  of  a  New¬ 
tonian  calorimeter;  however,  the  calorimeter  should  be  sized  so  that 
only  temperature  variation  along  one  axis  need  be  considered.  One  simple 
means  of  determining  the  severity  of  the  temperature  gradient  is  to 
assume  that  the  calorimeter  is  a  semi-infinite  solid  with  a  constant  heat 
rate.  Reference  11  presents  the  following  relation  for  the  temperature 
distributions  of  this  case  as  a  function  of  time  and  distance: 


Ta  (x,  6)  m  TAi  + 


Qa  ourc  e 
k  A 


~x2/a  a\Q 


—  x  erfc 


x 

2  V  ®A  0 


Now  for  any  given  time  (0)»  the  above  equation  can  be  solved  for  x  equal 
zero  and  x  equal  to  the  thickness  of  the  calorimeter,  and  the  resulting 
temperature  difference  can  be  analyzed.  It  must  be  remembered  that  the 
above  is  a  very  pessimistic  result  because,  for  the  application  in  this 
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report,  the  calorimeter  cannot  be  considered  as  a  semi-infinite  solid. 

If  the  calorimeter  cannot  be  considered  Newtonian,  there  are  several 
means  by  which  the  effect  of  the  temperature  gradient  may  be  compen¬ 
sated. 

1.  Construct  a  grid  network  similar  to  the  one  used  for  the  insulator. 

2.  Measure  the  surface  and  rear  face  temperature  of  the  calo¬ 
rimeter  and  assume  a  linear  relation  between  these  two  points. 

An  equation  of  the  following  type  is  obtained: 

where 

T(fl)t  is  the  surface  temperature  history 
T(0),  is  the  back  face  temperature  history 

3.  Embed  the  measuring  thermocouple  at  the  center  of  the  calo¬ 
rimeter  and  assume  the  resulting  temperature  history  is  the 
average  temperature  between  the  front  and  back  surface.  When 
embedding  a  thermocouple  within  the  calorimeter,  a  calculated 
estimate  should  be  made  of  the  error  caused  by  the  disturbance 
of  the  calorimeter  temperature  field  from  the  thermocouple 
insert. 


2.6  HEATING  OF  THE  BACK  SURFACE 

In  some  instances,  a  base  configuration  may  be  encountered  for 
which  it  is  necessary  to  determine  the  net  heating  rate  from  both  sides 
of  the  base,  such  as  a  protective  flame  shield  where  afterburning  of  fuel- 
rich  exhaust  products  may  occur  on  the  back  side.  A  calorimeter  con¬ 
figuration  as  shown  in  Fig.  3b  may  be  utilized;  however,  the  gas  film 
coefficients  cannot  be  determined  for  this  particular  design.  If  this  in¬ 
formation  is  desired,  two  calorimeters  as  shown  in  Fig.  3a,  inverted 
with  respect  to  each  other,  should  be  used. 


2.7  SURFACE  CONTAMINATION 

One  problem  encountered  with  the  use  of  the  calorimeter  is  the 
determination  of  any  change  in  surface  emissivity  which  may  be  caused 
by  rocket  exhaust  product  contamination.  One  method  which  has  been 
attempted  on  past  missile  base  heating  test  programs  was  to  use  a 
polished  calorimeter  (absorptivity  near  zero)  and  a  blackened  calorimeter 
(absorptivity  near  unity)  and  to  assume  that  the  difference  between  the 
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indicated  heat  rate  of  each  calorimeter  is  equal  to  the  source  radiant 
heat  rate.  Because  of  contamination  and  high  base  gas  temperatures, 
it  is  virtually  impossible  to  maintain  the  highly  reflective  calorimeter 
surface.  Fortunately,  the  indication  to  date  is  that  the  radiant  heat 
flux  is  generally  much  less  than  the  convective  heat  fluxes  for  most 
missile  bases.  Furthermore,  the  measuring  devices  proposed  in  this 
report  will  allow  the  radiant  heat  rate  to  be  measured  separately  from 
the  total  heat  rate,  and  this  will  assist  in  the  analysis  of  the  contamina¬ 
tion  problem  as  it  affects  the  slug  mass  calorimeter. 


3.0  ANALYSIS  OF  BLACK-BODY  CAVITY 


A  schematic  of  the  instrument  selected  to  measure  radiant  heat 
rates  is  shown  in  Fig.  4.  A  similar  device  is  reported  in  Ref.  12.  The 
principal  difference  is  that  the  instrument  used  in  Ref.  12  has  a  lens 
between  the  aperture  and  heat  sensing,  unit  in  conjunction  with  an  internal 
passage  through  which  air  is  bled  to  prevent  contamination  of  the  lens. 

In  addition,  the  instrument  in  Ref.  12  is  not  particularly  amenable  to 
mathematical  analysis  but  rather  relies  upon  interpolation  of  experi¬ 
mental  data  obtained  from  a  laboratory  calibration.  This  makes  it  neces¬ 
sary  to  maintain  the  external  environment  of  the  heat  sensor  similar  to 
that  for  the  calibration.  Although  the  use  of  a  lens  is  advantageous  in  that 
the  heat  sensing  unit  is  isolated  from  external  convection  currents,  the 
other  requirements  (air  jet  and  controlled  environment)  make  this  instru¬ 
ment  impractical  for  flight  test  missiles.  The  device  used  in  this  in¬ 
vestigation  has  an  open  aperture  through  which  the  radiant  heat  is  trans¬ 
mitted.  The  sizing  of  the  aperture  to  prevent  external  convective  currents 
from  strongly  influencing  the  cavity  temperature  history  is  a  problem  of 
an  extremely  complex  nature  (see  Fig.  5  for  a  pictorial  illustration  of  this 
effect);  therefore,  this  phenomenon  was  investigated  experimentally. 

3.1  MATHEMATICAL  MODEL 

In  analyzing  the  configuration  in  Fig.  4,  the  following  conditions  and 
assumptions  are  specified: 

1.  All  metal  parts  are  Newtonian. 

2.  The  thermocouple  leads  are  infinitely  long  and  perfectly  insulated. 

3.  Thermal  contact  resistances  at  interfaces  are  zero. 

4.  The  emissivity  of  the  cavity  is  equal  to  the  absorptivity. 

5.  Convective  heating  within  the  cavity  is  negligible. 
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6.  The  reflections  of  radiant  energy  within  the  cavity  are  diffused. 

7.  The  influence  of  the  metal  cap  on  the  temperature  gradient  in 
the  insulator  behind  the  heat  sensor  is  negligible. 

8.  All  of  the  incident  incoming  radiation  on  the  sides  of  the  aperture 
walls  is  specularly  reflected  (a  «<l)  and  no  energy  is  absorbed. 
This  further  implies  that  the  angle  factor  from  station  1  to 
station  2  (see  Fig.  6)  is  unity. 

9.  The  absorptivity  for  the  surface  of  the  heat  sensor  and  the  back 
surface  of  the  metal  cap  is  near  unity. 

10.  The  insulating  washer  which  separates  the  metal  cap  and  heat 
sensor  has  negligible  effect  on  angle  factors  with  respect  to  the 
heat  sensor. 

11.  The  temperature -time  history  for  the  heat  sensor  may  be  ex¬ 
pressed  by  the  following  empirical  relation: 

TD(0)  -  Ti  +  BMl  -  (6) 

where 

B'  -  experimental  constant,  °R 
Z'  -  experimental  constant,  per  hr 

12.  The  temperature -time  history  of  the  metal  cap  may  be  expressed 
by  the  following  empirical  relation: 

T(0)  -  Ti  +  CMl  -  fTE'9)  (7) 

where 

C'  -  experimental  constant,  °R 
E'  =  experimental  constant,  per  hr 

For  any  given  time,  the  following  heat  balance  may  be  written  about  the 
heat  sensor  (refer  to  Fig.  6): 

Qd  (1  -*  5)  +  Qd(4-*5)  =  QdI  +  QD2  +  QD3  +  QD4  (8) 

where 

Qd  (1  -*5)  =  radiation  through  the  aperture  to  the  heat 
sensor 

Qd  (4  -»5)  =  re -radiation  from  the  base  of  the  metal 
cap  to  the  heat  sensor 
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Qdi  =  change  in  internal  energy  for  the  heat  sensor 
Qd2  =  conduction  losses  through  the  thermocouple  leads 
Qd3  =  conduction  losses  into  the  insulation 
Qd4  =  re -radiation  from  the  heat  sensor 

Furthermore, 

Qd  ( 1  -*  5 )  =  «0  AD  F,_s  a  Ts  (0)4  (8a) 

where 

*•  -  ..(.--■fcrri7rl>  +  (1 ' '■>' u/s  - •*•>] (Ref'  13) 

Fi-.  -  1 

Ts  ( 0 )  =  apparent  radiant  source  black-body  temperature 
Qd(4-*5)  =  e4A4F4_5  a  [c '  (1  -  )  +  T|]  (8b) 


where 

F4_,  =  1 

Qd4  -  «dADFD  a[B'(l  -e~Z'd  )  +  Ti] 


where 


Fq  =  sin  ijj 


Qdi 


n*D  cpp 


B'Z' 


-z'd 

e 


QD2 


-kc  ttt 


2 

c 


dTp  (d) 
dxC 


*c 


0 


QD3 


-kB  n  V  +£  1 


d  Tp  (0)1 


(8c) 

(Ref.  14) 

(8d) 

(8e) 

(8  f) 
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3 2  SOLUTION  OF  THERMOCOUPLE  LOSSES 


Equation  (8e)  can  be  solved  by  using  Eq.  (3),  and  the  resulting 
equation  is 

QD2  .  ^crc*  2B'Z'£1/a  t~Z'd  +  -l'[l  T-* - —  - - —  +  JLL 

V  ac/ff  /3‘4  e% 

-  +  2Z'i8,/a  -  -1  Z'2  /3V*  +  -JL_Z'*  /3Va  -  ■  •  •] 


+  2Z'  el/l _ 1  Z'2  efi  +  z'2  0*4  -  . . 

3  15 


for  /3  >  0  ,  Z '  >1,  0  <  1 


3.3  SOLUTION  OF  INSULATOR  LOSSES 


If  the  insulation  can  be  considered  as  a  one -dimensional  semi¬ 
infinite  slab  with  a  uniform  initial  temperature,  Eq.  (8f)  may  also  be 
solved  by  Eq.  (3),  and  the  resulting  equation  is 

qD3  .  r2B-z- e‘4 ,-«••  *  B'll-l/<9~gl1  -  ♦  4 

V  «c/»T  1  P'*  #4 

-  B'  e”2*  (0-0)  T - 1  +  2Z'  fi1/a - i-  Z'a  Rl/i  +  — Z'2  jS*^  -  .  .  .]  (10) 

L  /3/j  3  15  J 

+  B'  T _ L_  +  2Z'0l/j  -  i  Z'1  5^  +  Z'*  5 -  •  •  •] 

L  .y.  3  15  J 


for  >3  >  0  ,  Z'>1,  0  <  1 

To  determine  the  amount  of  insulation  which  is  required  to  satisfy 
Eq.  (10),  the  following  equations  which  are  presented  in  Ref.  11  may 
be  used. 

Boundary  conditions  for  the  heat  sensor  are 
Tb  -  Ti  s  0  =  0  ,  x  'B  >  0 

TB  =  Ti  +  B'(l  -e“Z'*)  «  x'B  =0j  6>  o 

Tb  =  Ti  e  x'B  =  ~,0*  0 
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The  solution  is 


^2 *  Q 

Tb(<9,x#b)  =  Ti  +  Bl  erfc  XB  +  -^-2 - 


+  eX  B  ‘Vz'/as  erfc  XfllJ 


|e“X  B  W  Z'/“B  erfc  Xbo 


(ID 


where 

XB 

Xbo 

Xbi 


2  V  aB  0 

XB  -  i  VZ'0 

xB  +  iVz'0 


The  same  equation  may  be  used  for  the  external  casing.  Tables  for  the 
solution  of  probability  integrals  for  complex  argument  are  presented  in 
Ref.  15.  Now  using  Eq.  (11),  a  temperature  history  for  the  heat  sensor 
and  external  casing  may  be  assumed,  and  the  depth  of  penetration  for  an 
allowable  temperature  rise  for  any  given  time  determined. 


3.4  HIGH  RADIANT  HEAT  RATE  VALUES 

For  relatively  high  radiant  heat  rates,  it  is  not  desirable  to  isolate 
the  heat  sensor  from  the  external  casing  because  the  large  temperature 
gradient  in  the  insulator  will  void  the  use  of  constant  thermal  properties 
when  the  heat  losses  into  the  insulator  are  calculated.  Although  the  in¬ 
sulator  losses  may  be  calculated  by  determining  the  thermal  properties 
as  a  function  of  the  temperature  level  at  given  depths  in  the  insulation,  it 
is  felt  that  the  instrument  shown  in  Fig.  7  will  simplify  this  procedure. 

The  insulation  losses  for  this  instrument  are  determined  by  using  a  two- 
dimensional  grid-network  (as  previously  discussed  in  section  2.  3)  with 
the  thermal  diffusivity  of  the  insulator  determined  at  the  temperature  level 
of  the  heat  sensor  or  by  using  the  average  temperature  between  the  heat 
sensor  and  external  casing,  whichever  is  more  representative  of  the  tem¬ 
perature  level  in  the  insulation.  It  may  not  always  be  possible  to  assume 
that  the  external  casing  is  Newtonian,  but  this  situation  can  be  handled  in 
the  same  manner  as  previously  discussed  for  non-Newtonian  slug  mass 
calorimeters. 


3.5  LOW  RADIANT  HEAT  RATE  VALUES 

Since  the  ratio  of  the  surface  area  of  the  heat  sensor  to  the  aperture 
area  is  relatively  large  for  the  black-body  cavities  considered,  radiant 
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heat  fluxes  less  than  approximately  3  Btu/ft^-sec  will  be  difficult  to 
determine  accurately  because  of  the  small  temperature  rise  of  the  heat 
sensor.  Fortunately,  for  the  purpose  of  base  heating  tests,  heat  fluxes 
of  this  magnitude  are  insignificant  and  can  generally  be  ignored.  In 
other  type  testing  applications  it  may  be  desirable  to  measure  these 
values.  The  black-body  configuration  shown  in  Fig.  8  is  analyzed  as  a 
means  of  measuring  the  lesser  rates  of  radiant  heat  fluxes.  The  prin¬ 
cipal  advantage  of  this  instrument  is  that  the  air-gap  between  the  heat 
sensor  and  backing  plate  may  be  evacuated,  and  the  conduction  losses 
from  the  back  of  the  heat  sensor  are  greatly  reduced.  The  conditions 
and  assumptions  for  this  design  are  the  same  as  previously  stated  for  the 
other  black-body  cavity  with  the  following  additions: 

1.  The  pressure  in  the  evacuated  region  between  the  heat  sensor  and 
the  backing  plate  is  maintained  in  a  rarefied  regime  so  that  ther¬ 
mal  conductance  in  this  region  is  negligible. 

2.  The  angle  factor  between  the  heat  sensor  and  backing  plate  may 
be  determined  by  considering  the  two  bodies  as  infinite  parallel 
planes  each  at  uniform  temperature  at  any  time  (0) 

3.  The  temperature  history  for  the  heat  sensor  and  backing  plate 
may  be  expressed  by  the  same  type  exponential  equation  as  pre¬ 
viously  used  for  the  heat  sensor. 

Now  the  following  heat  balance  may  be  written  about  the  heat  sensor  (re¬ 
fer  to  Fig.  6): 

Qd  (1  -»  5)  +  Qd  (4  -*  5)  =  Qdx  +  QDj  +  Qd4  +  Qd  (5  %  6)  ^2) 

where 

Qd  (5  $  6)  =  radiation  interchange  between  the  heat 
sensor  and  the  backing  plate 

also 

Qd(s;6)  .  1/,;i;,;;-|[TDW-I,W]  (13) 

where 

To  (0)  =  temperature  heating  for  the  heat  sensor 

T«  (0)  =  temperature  history  for  the  backing  plate 

The  other  terms  in  Eq.  (12)  are  the  same  as  those  used  in  Eq.  (8).  It 
should  be  noted  that,  for  low  emissivity  values  for  the  back  of  the  heat 
sensor  and  the  face  of  the  backing  plate,  the  corresponding  temperature 
rise  for  the  heat  sensor  at  a  given  source  heat  rate  is  much  greater. 
Therefore,  polished  surfaces  are  desirable  because  the  equation  which 
will  express  the  temperature  history  of  the  heat  sensor  can  be  more 
accurately  defined. 
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3.6  MEASUREMENTS  IN  SPECIFIC  AREAS 

One  advantage  to  the  use  of  an  aperture  is  that  it  may  be  sized  to 
measure  radiant  heat  rates  from  specific  areas  in  the  base  region  if 
the  side  walls  of  the  aperture  are  made  of  a  non- reflecting  material 
(absorptivity  is  unity)  so  that  only  the  direct  radiation  will  enter  the 
cavity. 

The  angle  factor  for  direct  radiation  is  (refer  to  Pig.  6) 

Fj_j  =  1+2  (-A-)*  -  2  -i-  V(x/D)a  +  1  (Ref.  16)  (14) 

Therefore,  by  properly  sizing  the  x/D  ratio  for  the  aperture,  the  field 
of  view  for  the  cavity  can  be  controlled.  In  blackening  the  aperture 
walls,  the  re -radiation  from  the  walls  of  the  aperture  to  the  heat  sen¬ 
sor  may  become  significant,  and  this  may  be  accounted  for  as  follows 
(refer  to  Fig.  6) 

QD  (3-5)  =  *  A,  Fj_,  a  [Ti  +  C'(l  -  e“E  *)]  (15) 

where 

F»-s  =  lA  [Vfx/D)2  +  Y  -  x/D]  (Ref.  14) 


4.0  EXPERIMENTAL  PROGRAM 


4.1  CALIBRATION  FURNACE 


To  determine  the  effectiveness  of  the  general  mathematical  pro¬ 
cedures  as  used  in  this  report,  a  series  of  bench  calibrations  was  con¬ 
ducted  with  a  black-body  radiant  energy  source.  The  particular  design 
used  for  the  black-body  was  adapted  from  the  one  described  in  Ref.  17 
and  illustrated  in  Fig.  9.  The  inner  liner  which  simulates  the  black-body 
is  graphite  (f  =  0.95),  and  by  using  the  following  equation 


where 


£  gr  (1  —  s/S )  +  s/S 


1  +  (1  - 


<gr)(s/S  -  s/S0 ) J 


£gr  =  emissivity  of  the  graphite 


(16) 


the  emissivity  of  the  liner  aperture  is  calculated  to  be  greater  than  0.  99. 
A  description  of  the  furnace  used  for  the  calibrations  is  as  follows: 

Heating  elements  -  molybdenum 
Insulation  -  firebrick 
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Maximum  operating  temperature  -  3560°R 

Maximum  temperature  variation  -  ±4  deg 
within  the  center  section 


Inert  atmosphere 
Helium  flow  rate 


-  helium 

-  approximately 
80  scfh 


The  temperature  of  the  furnace  was  automatically  controlled  by  a  photo¬ 
electric  cell  which  was  sighted  on  the  furnace  wall.  The  temperature  of 
the  graphite  liner  was  measured  with  a  planinum-platinum- rhodium 
thermocouple  monitored  on  a  sixteen-point  recorder  and  also  with  an 
optical  pyrometer.  The  calibration  of  the  pyrometer  was  compared  to 
a  National  Bureau  of  Standards  calibrated  tungsten  filament.  All  other 
temperatures  were  recorded  on  an  eighteen- channel  oscillograph.  The 
uncertainty  of  the  source  heat  rate  calculated  from  the  furnace  tempera¬ 
ture  measurements  is  estimated  to  be  five  percent. 


To  determine  the  magnitude  of  the  thermal  effect  of  external  forced 
convection  on  the  black-body  cavity,  an  inconel  tube  was  inserted  into  the 
graphite  liner  to  permit  flowing  of  helium  gas  over  the  surface  of  the 
instrument  (Fig.  10).  The  intent  of  this  investigation  was  not  to  bracket 
the  phenomenon  completely  but  only  to  determine  the  strength  of  this 
effect  for  external  velocities  which  were  felt  to  be  representative  of  the 
upper  values  encountered  on  missile  bases  and  for  cavity  apertures  not 
exceeding  0.  15  in.  Because  of  uncertainty  in  the  exact  temperature  rise 
of  the  helium  gas  as  it  passed  through  the  inconel  tube,  the  exit  velocity 
of  the  gas  from  the  tube  was  conservatively  estimated,  based  on  the  volume 
flow  and  pressure  measurement  of  the  gas,  to  be  1500  ft/sec.  The  leakage 
rate  of  the  graphite  liner  and  furnace  was  great  enough  so  that  for  steady- 
state  conditions  the  liner  was  not  pressurized. 


4.2  SLUG  MASS  CALORIMETER 

The  slug  mass  calorimeter  used  for  the  program  is  illustrated  in 
Fig.  11.  Iron-constantan  thermocouples  (28  gage)  were  spot  welded  to 
the  slug  mass  and  to  the  side  and  back  of  the  metal  casing.  The  insulating 
material  was  a  polytrifluorochloroethylene  plastic.  The  calorimeters  were 
stainless  steel  type  347,  and  each  was  weighed  and  measured  prior  to 
installation.  The  metal  casings  for  each  instrument  had  a  45 -deg  lip  to 
ensure  that  the  surface  of  the  calorimeter  was  flush  with  the  inside  of  the 
graphite  liner  and  had  a  rod  attached  to  the  rear  of  the  casing  to  enable 
insertion  into  the  liner  aperture.  Information  was  available  for  the 
emissivity  of  the  calorimeter  («  =  0.39  for  T  <  860=11,  Ref.  18).  The 
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surface  emissivity  of  the  insulator  was  unknown,  therefore,  the  insu¬ 
lator  was  coated  with  a  rough  layer  of  lamp  black  (e  =  0.8,  Ref.  18). 


4.3  BLACK-BODY  CAVITY 

A  schematic  of  the  black-body  cavity  which  was  calibrated  is  shown 
in  Fig.  12.  The  entire  internal  surface  of  the  cavity  was  coated  with  a 
commercial  optical  black  paint,  while  the  side  walls  of  the  aperture 
were  polished.  Iron-constantan  thermocouples  (28  gage)  were  spot  welded 
to  the  copper  heat  sensor  and  the  stainless  steel  surface  cap.  The  insu¬ 
lating  material  was  polytrifluorochloroethylene  plastic.  Each  of  the  copper 
cones  was  weighed  and  measured  prior  to  installation.  The  aperture  sizes 
which  were  used  for  this  program  were: 

Aperture  Diameter,  in.  Aperture  D/x 

0.05  1.67 

0.15  5.00 

All  of  the  thermal  properties  and  physical  dimensions  used  in  this 
report  are  presented  in  Tables  1  and  2  and  Figs.  13  through  19. 


4.4  PROCEDURE 

1.  The  desired  furnace  temperature  was  set  with  the  automatic 
control  system,  and  after  obtaining  the  pre-set  temperatures  the 
furnace  was  allowed  to  soak  for  an  additional  hour  to  ensure  tem¬ 
perature  uniformity. 

2.  After  the  soaking  period  had  been  completed,  the  optical  pyrom¬ 
eter  was  sighted  through  the  front  of  the  furnace  in  order  to 
determine  accurately  the  graphite  liner  temperature. 

3.  The  flow  of  the  inerting  atmosphere  was  reduced  to  less  than 
10  scfh. 

4.  The  enclosure  over  the  front  of  the  aperture  in  the  furnace  door 
was  removed  and  the  instrument  inserted  into  the  graphite  liner 
aperture.  The  time  required  to  seat  the  instrument  was  on  the 
order  of  one  second.  When  the  instrument  was  seated,  a  blip  on 
the  oscillograph  trace  indicated  time  zero. 

5.  When  determining  the  thermal  effect  of  external  forced  convection 
on  the  black-body  cavity,  the  same  procedure  as  described  above 
was  followed  with  the  addition  that  when  the  instrument  was  seated 
in  the  graphite  liner  an  off- on  hand  valve  in  the  external  flow  simu¬ 
lation  tube  was  opened. 
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5.0  RESULTS  AND  DISCUSSION 
5.1  SLUG  MASS  CALORIMETER 

Two  typical  experimental  temperature -time  histories  for  the  calorim¬ 
eter  along  with  the  respective  boundary  temperature  histories  are  illus¬ 
trated  in  Figs.  13a  and  b  with  the  corresponding  exponential  equation 
presented  with  each  curve.  In  fitting  the  exponential  equation  to  the  ex¬ 
perimental  temperature -time  curve,  the  deviation  was  never  more  than  a 
few  degrees;  therefore,  only  the  experimental  curve  was  presented.  The 
respective  corrected  heat  rate-time  histories  are  presented  in  Figs.  14a 
and  b.  Change  in  internal  energy,  re -radiation,  and  thermocouple  losses 
were  calculated  as  outlined  in  section  2.1.  An  emissivity  value  of  0.  4 
was  used  in  determining  the  radiation  losses  from  the  calorimeter  surface. 
It  was  first  considered  that  the  calorimeter  emissivity  might  change  sig¬ 
nificantly  when  exposed  to  the  furnace  temperature;  however,  it  was  found 
that  no  significant  change  occurred  in  the  temperature-time  history  of  the 
calorimeter  when  it  was  used  twice  to  measure  the  same  source  heat  rate. 
The  fact  that  there  was  no  apparent  change  in  the  calorimeter  surface 
emissivity  is  attributed  to  the  inerting  helium  atmosphere  and  the  short 
exposure  time.  Also,  each  calorimeter  was  used  only  a  limited  number 
of  times.  Equation  (4),'  which  is  used  to  determine  the  heat  loss  into  the 
thermocouple  leads,  was  programmed  on  a  7070  EBM  digital  computer 
and  the  series  which  is  part  of  this  equation  was  expanded  to  include  twenty 
terms.  In  determining  the  losses  and  gains  from  and  into  the  insulation, 
a  grid-network  45  units  in  depth  and  50  units  in  width  was  used.  The  heat 
balance  equations  for  the  grid-network  were  also  programmed  on  the 
7070  IBM  digital  computer,  and  the  temperature  field  in  the  insulation  for 
specified  time  periods  was  obtained  as  the  end  result.  The  grid  network 
was  solved  in  the  same  manner  as  previously  discussed  except  QBTotal  for 
this  instance  was: 


QBTotml 


£ B  St<r( Tf4  -  T;4  ) 


where 


Tf  =  graphite  liner  temperature,  °R 

Tj  =  surface  grid  temperature  after  time  A0,  °R 

The  temperature  of  the  surface  grids  was  calculated  by  iterating  the 
individual  grid  heat  balance  until  the  temperature  determined  from  the 
incoming  radiant  heat  rate  was  equal  to  the  corresponding  temperature 
for  the  heat  balance  at  the  particular  time  increment.  The  tabulated  tem¬ 
peratures  from  the  digital  computer  were  plotted  graphically  versus  posi¬ 
tion  with  respect  to  the  calorimeter.  Figures  15a  and  b  show  two  typical 
temperature  distributions  in  the  insulator  adjacent  to  the  side  and  rear. 
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respectively,  of  the  calorimeter.  The  net  exchange  of  heat  between  the 
calorimeter  and  insulator  was  determined  by  using  a  planimeter  to  meas¬ 
ure  the  area  between  the  calorimeter  temperature  curve  and  the  tem¬ 
perature  distribution  in  the  insulator  next  to  the  calorimeter.  In  retro¬ 
spect,  because  of  the  upper  temperature  limit  of  flurocarbon  plastics,  a 
ceramic-type  material  would  have  served  more  favorably  as  the  insu¬ 
lating  material;  however,  in  consideration  of  the  relatively  short  time 
periods  involved,  no  significant  adverse  effects  were  noticeable.  In 
addition,  the  thermal  characteristics  of  this  material  made  it  possible  to 
approximate  constant  values  for  the  thermal  conductivity  and  specific 
heat  at  the  average  temperature  of  the  calorimeter. 

With  reference  to  Figs.  14a  and  b,  the  curvature  of  the  calculated 
source  heat  rate  is  the  effect  of  using  constant  thermal  properties  for 
the  insulator  in  the  calculation  of  the  grid-network.  This  curvature  will 
become  more  pronounced  as  the  source  heat  rate  increases,  and,  in 
actual  testing  applications  the  change  of  thermal  properties  with  tempera¬ 
ture  must  be  considered.  The  important  factor,  therefore,  in  selecting 
an  insulating  material  for  use  with  the  calorimeter  is  that  it  have  known 
thermal  properties  as  a  function  of  temperature.  Reference  19  presents 
the  thermo-physical  properties  as  a  function  of  temperature  for  numerous 
insulating  materials.  Figures  14a  and  b  also  illustrate  that  for  even  a 
good  insulation  material  the  conduction  losses  into  the  insulator  can  be 
quite  significant.  In  fact,  the  calorimeters  used  in  base  heating  tests  will, 
in  all  probability,  be  thicker  and  exposed  to  heating  rates  greater  than  the 
values  encountered.  These  factors  will  further  increase  the  error  in 
assuming  that  the  source  heat  rate  is  equal  to  the  change  in  internal  energy 
of  the  calorimeter. 

The  results  of  the  mathematical  computations  illustrate  that  total 
source  heat  rates  can  be  calculated  with  reasonable  accuracies  if  a  com¬ 
plete  heat  balance  is  written  about  the  calorimeter.  Also,  the  accuracy 
in  determining  .source  heat  rates  should  improve  as  the  source  heat  rate 
increases  because  a  small  error  in  the  measurement  of  the  temperature 
or  in  curve  fitting  will  have  a  less  pronounced  effect  on  the  results. 

5.2  BLACK-BODY  CAYITY 

Two  typical  temperature -time  curves  for  the  black-body  cavity  with¬ 
out  external  flow  are  presented  in  Figs.  16a  and  b.  The  empirical  equation 
used  for  the  purpose  of  calculations  is  presented  with  each  curve.  The 
respective  heat  rate-time  histories  are  shown  in  Figs.  17a  and  b.  The 
change  in  internal  energy  and  re-radiation  gains  and  losses  were  calculated 
as  outlined  in  section  3.  1.  An  emissivity  value  of  unity  was  used  in  deter¬ 
mining  the  radiation  from  the  heat  sensor  and  surface  cap.  The  radiation 
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interchange  with  the  sides  of  the  aperture  wall  was  assumed  to  be  neg¬ 
ligible.  The  thermocouple  and  insulator  losses  for  the  experimental 
curve  in  Fig.  16b  were  calculated  from  Eqs.  (9)  and  (10),  respectively 
In  determining  the  thermocouple  and  insulator  losses  for  the  experi¬ 
mental  curve  in  Fig.  16a  the  following  equation  was  used: 

Q  =  2  lc  A  c  6V*  (17) 

\Jn  a 

where 

C  is  the  temperature -time  slope 

of  the  experimental  curve,  deg/hr 

Equation  (17)  is  the  solution  of  the  one -dimensional  heat  conduction 
equation  of  a  semi- infinite  slab  with  a  uniform  initial  temperature 
where  the  surface  temperature  history  is  a  linear  relation  with  time. 

The  initial  thermal  lag  of  the  heat  sensing  unit  (Figs.  16a  and  b)  is 
attributed  to  the  magnitude  of  the  heat  sensor  surface  area  compared 
to  the  magnitude  of  the  heat  rate  encountered  and  to  the  location  of  the 
thermocouple  leads  at  the  apex  of  the  cone.  Therefore,  in  calculating 
the  change  in  internal  energy,  insulator  losses,  and  thermocouple 
losses,  time  zero  was  considered  as  the  point  in  the  initial  portion  of 
the  temperature -time  history  where  an  inflection  occurred  in  the  slope 
of  the  curve.  In  Figs.  16a  and  b  the  inflection  was  considered  to  be  at 
3.  75  and  2.  5  sec,  respectively.  This  inflection  point  may  be  difficult 
to  determine  accurately,  particularly  for  low  magnitudes  of  radiant 
heat  rates.  The  source  heat  rate  was  calculated  for  the  temperature 
curve  in  Fig.  16b  assuming  that  the  inflection  point  occurred  at  1.  25  sec 
rather  than  2.  5  sec.  The  temperature -time  relation  then  becomes 

Ta(0)  -  205  (1  -  e“780)  +  541.5 

and  the  heat  rate-time  history  using  this  expression  is  shown  in  Fig.  17c. 
A  comparison  of  Figs.  17b  and  c  shows  that  determination  of  the  zero 
time,  within  reason,  which  corresponds  to  the  inflection  point,  does  not 
have  an  appreciable  effect  on  the  calculated  source  heat  rate. 

Two  typical  temperature -time  curves  for  the  black-body  cavity  with 
external  flow  are  illustrated  in  Figs.  18a  and  b.  The  exponential  equation 
used  for  the  calculations  is  presented  with  the  heat  sensor  curve.  The 
re-radiated  heat  rate  from  the  surface  cap  to  the  heat  sensor  was  calcu¬ 
lated  using  the  measured  temperatures  because  of  the  difficulty  of  fitting 
an  accurate  exponential  equation  to  this  curve.  The  respective  calculated 
heat  rate-time  histories  for  the  temperature  curves  in  Figs.  18a  and  b 
are  presented  in  Figs.  19a  and  b.  The  results  in  Figs.  19a  and  b  indicate 
that,  initially,  the  black-body  cavity  was  subjected  to  very  high  heating 
rates  and  then  approached  equilibrium  at  the  radiant  source  heat  rate 
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received  from  the  graphite  liner.  Because  the  calculated  source  heat 
rate  in  Figs.  19a  and  b  are  still  in  a  transient  state  at  the  end  of  the 
measured  time  period,  an  equation  of  the  following  form  was  fitted  to 
the  calculated  source  heat  rate 

Q(d)  =  A  +  B  e~zd 

and  d  was  taken  to  infinite.  The  equilibrium  calculated  source  heat 
rate  for  Figs.  19a  and  b  was,  respectively,  zero  and  600  x  10_6Btu/sec. 
The  magnitude  of  the  error  for  the  calculated  equilibrium  versus  the 
actual  source  heat  rate  for  Figs.  19a  and  b  is  not  too  disconcerting  when 
considering  the  accuracy  necessary  to  obtain  the  calculated  equilibrium 
source  heat  rate  after  experiencing  such  a  high  initial  value.  The  high 
initial  heat  rate  is  not  attributed  to  a  mass  interchange  process  as  illus¬ 
trated  in  Fig.  5  because  it  is  felt  that  if  this  phenomenon  were  predom¬ 
inating,  the  heat  sensor  would  have  obtained  much  higher  temperatures, 
and  the  calculated  equilibrium  source  heat  rate  would  have  had  a  value 
much  greater  than  the  actual  source  radiant  heat  rate.  Rather,  the  high 
initial  heating  rates  in  Figs.  19a  and  b  are  attributed  to  a  large  pressure 
surge  of  helium  gas  at  near  furnace  t'emperature  into  the  cavity  when  the 
valve  in  the  external  flow  simulation  tube  was  opened.  This  is  further 
evidenced  by  the  fact  that  the  initial  temperature  rise  for  both  instru¬ 
ments  is  independent  of  the  fact  that  the  cavity  aperture  diameters  were 
0.  05  and  0.  15  in.  in  diameter.  Fortunately,  in  testing  at  a  given  simu¬ 
lated  altitude  condition,  the  source  radiant  heat  rate  should  remain 
relatively  constant.  Therefore,  if  a  pressure  surge  involving  hot  gases 
is  encountered  during  the  start  transient  of  a  rocket  test,  the  black-body 
cavity  calculated  source  heat  rate  can  be  solved  at  greater  times  until 
a  constant  valve  is  obtained.  Furthermore,  for  purposes  of  actual  flight 
tests,  in  general,  missile  base  pressures  initially  decrease  slowly  and 
then  tend  to  remain  relatively  constant  with  time. 

One  method  of  improving  the  current  black-body  design,  as  evidenced 
by  the  experimental  data,  would  be  to  replace  the  metal  surface  cap  with 
a  cap  made  of  an  insulating  material,  such  as  a  ceramic  and  to  use  a 
polished  metallic  insert  as  the  cavity  aperture.  This  change  would  serve 
a  two-fold  purpose  in  that  it  would  eliminate  an  extraneous  source  of  heat 
to  the  heat  sensing  unit  and  also  would  limit  internal  convective  currents 
caused  by  temperature  variations  within  the  cavity. 

In  all,  an  aperture  diameter  of  0.  15  in.  may  be  used  in  conjunction 
with  the  black-body  cavity  and  is  more  desirable  than  the  0.  05-in. - 
aperture  because,  for  the  same  cavity  internal  surface  area,  the  resulting 
temperature-time  history  for  the  same  source  heat  rate  will  be  greater 
with  the  0.  15 -in.  aperture.  Thus,  the  empirical  equation  used  to  define 
the  temperature  history  can  be  more  accurately  defined,  and  also  the 
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absolute  effect  of  the  instrument  read-out  error  will  be  lessened.  One¬ 
dimensional  heat  conduction  equations  may  be  used  to  determine  insu¬ 
lator  losses  providing  sufficient  insulation  is  used  based  on  the  environ¬ 
mental  temperature  histories  and  the  exposure  time  period.  However, 
if  so  desired,  it  should  be  possible  to  improve  the  accuracy  of  the  results 
by  calculating  the  insulator  losses  by  the  method  of  finite  differences  and 
accounting  for  the  emissivity  of  the  sides  of  the  aperture  walls. 


6.0  CONCLUDING  REMARKS 


Some  of  the  methods  of  measuring  heat  rates  in  the  base  areas  of 
rocket-powered  missiles  and  the  inherent  errors  in  these  methods  are 
summarized  as  follows: 

1.  The  use  of  a  polished  and  blackened  slug  mass  calorimeter  to 
determine  convective  and 'radiant  heat  rates  can  result  in  con¬ 
siderable  error  because  of  surface  contamination  on  the  polished 
surface  caused  by  the  rocket  exhaust  gases. 

2.  The  practice  of  embedding  a  thin  disc  copper  calorimeter  in  a 
steel  missile  base  can  cause  considerable  error  because  of  the 
large  change  in  the  base  temperature -time  history  as  compared 
to  the  temperature-time  history  of  the  segment  it  replaces. 

3.  When  analyzing  a  slug  mass  calorimeter,  the  assumption  that 
the  incoming  heat  rate  is  equal  to  the  change  in  internal  energy 
may  be  greatly  in  error  because  of  the  resultant  gains  and/or 
losses  of  heat  through  the  surrounding  insulation. 

4.  The  practice  of  calibrating  a  device  which  will  measure  convective 
heat  rates  in  a  laboratory  with  a  standard  black-body  radiant 
energy  source  and  then  comparing  experimental  results  with  the 
calibration  to  determine  heat  rates  requires  that  the  external 
environment  of  the  calorimeter  be  artifically  controlled  or  isolated 
in  order  to  duplicate  the  conditions  during  the  calibration.  Either 
of  these  methods  will  cause  a  considerable  change  in  the  actual 
base  thermal  characteristics  with  a  resultant  error  because  of  the 
effects  on  the  base  gas  film  coefficient. 

In  this  report  it  has  been  determined  that  an  embedded  slug  mass 
calorimeter  and  a  black-body  cavity  may  be  used  to  obtain  reasonably  accu¬ 
rate  heat  transfer  measurements  in  the  base  areas  of  rocket-powered 
missiles  if  an  analytic  heat  balance  is  written  about  the  heat  sensing  unit. 
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When  using  these  devices  to  measure  heating  rates,  the  following  practices 
should  be  observed: 

1.  Preliminary  analysis  to  properly  size  and  design  the  measuring 
devices.  (One  of  the  prime  considerations  in  the  design  should  be 
to  facilitate  the  mathematical  analysis. ) 

2.  Measurement  of  each  instrument  to  determine  the  heat  transfer 
areas  required  for  the  mathematical  computations. 

3.  Individual  weight  determination  of  the  calorimeter  and  cavity  heat 
sensing  unit. 

4.  Extreme  care  in  instrumentation  and  data  reduction  procedures. 
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TABLE  1 

THERMAL  PROPERTIES  USED  FOR  THE  EXPERIMENTAL  PROGRAM 


Material 

Density, 

lbm/ft3 

Conductivity, 

Btu/ ft2-hr-deg/ft 

Specific  Heat, 
Btu/lbm-deg 

Polytrifluorochloroethylene 

Plastic 

135 

0.  035  -  0.  07 
at  (550  -  780°R) 

0.  22  -  0.  27 
at(550  -  780°R) 

Iron 

492 

36 

0.  106 

Constantan 

555 

13 

0.  102 

Stainless  Steel 

Type  347 

0.  13 

Copper 

0.  094 

TABLE  2 

PHYSICAL  CONSTANTS  USED  FOR  THE  BLACK-BODY  CAVITY 


Aperature, 

in. 


Area  (Ai), 
ft2 


Area  (A4), 
ft2 


Surface  Area  (Ad), 
ft* 


Diameter  (Di), 
ft 


Heat  Sensor, 
Mass,  grams 


0.  05 

1.  363  x  10-5 

3.271  x 

IO-4 

1.  6  x 

IO-3 

2.  083  x  10"  2 

0.  15 

1.  2266  x  10-4 

2. 1805  x 

IQ"4 

1.  6  x 

IQ'3 

2. 083  x  10"  2 

0.  597 
0.  582 


AEDC-T  DR-63-93 


AEDC-TDR.63-93 


QA,  Total 


Calorimeter 


Side  View 


Fig.  1  Slug  Mass  Calorimeter 
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AEDC-TDR- 63-93 


a.  Surface  Heat  Transfer  Measurements 


b.  Surface  and  Back  Heat  Transfer  Measurements 
Fig.  3  Slug  Mass  Calorimeters 
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AEDC-TDR-63-93 
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Thermocouple  Leads 


Fig.  4  Black-Body  Cavity 
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Fig.  5 


External 


Flow  Fi. 


•3 


AEDC-TDR.63-93 
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Vacuv 

Tube 


* 


ermocouple 
Leads 


Fig.  8  Black-Body  Cavity  (Low  Hoat  Ratos) 
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Fig.  9  Experimental  Black-Body  Radiation  Source 


Fig.  10  External  Forced  Convection  Simulation  Tube 
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AEDC-TDR-63*93 


Fig.  11  Experimental  Slug  Mass  Calorimeter 
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Thermocouple  Lead 


All  Dimensions  in  Inches 


Fig.  12  Experimental  Black-Body  Cavity 
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AEDC-TDR.63-93 


Calorimeter  Mass  «  0.811  grams 
Furnace  Temperature  -  1960°R 


900 


800 


+i 

a 

u 

& 

§ 


700 


600  m 


500  5 


10  20 
Time,  0,  sec 

a-  Qsouree  =  7.02  Btu/ft2 ■  sec 


Fig.  13  Calorimeter  Temperature-Time  Histories 
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AEDC-TOR-63-93 


cf 

■» 

0) 

a 

•p 


Calorimeter  Mass  ■  0.864  grams 
Furnace  Temperature  <-  2310°R 


Time,  9,  sec 


b.  Q.ourc.  =  13.55  Btu/ft2-ioc 
Fig.  13  Concluded 
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Heat  Rate  x  10  ,  Btu/sec 


AEDC-TDR-63-93 


-4  2 

Surface  Area  of  Calorimeter  -  7.466  x  10  ft 
Kg  -  0.05  Btu/ft2-hr-deg/ft 
CpB  -  0.25  Btu/lbffl-deg 
— — —  Calculated  Source  Heat  Rate 


Time,  0,  sec 


Qsource  -  7.02  Btu/ft2  •sec 
Fig.  14  Colorimeter  Heot  Rote-Time  Histories 
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AEDC-T  DR-63-93 


-4  2 

Surface  Area  of  Calorimeter  -  7 . 466  x  10  ft 
Kg  -  0.07  Btu/ft^-hr-deg/f t 


C  _  -  0.27  Btu/lb  -deg 
pc  m 

Calculated  Source  Heat  Rate 


b.  Qsource  -  13.55  Btu/ft2-S»C 
Fig.  14  Concluded 
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Temperature 


AEDC*TDR>63*93 


Distance,  x (0.574)  x  10“^,  ft 
a.  Side  Temperature  Distribution 

Fig.  15  Typical  Temperature  Distributions  in  the  Insulator  Adjacent  to  the  Calorimeter 
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Temperature 


AEDC*TDR>63-93 
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b.  Bock  Temperature  Distribution 
Fig.  15  Concluded 
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AEDC-TDR-63-93 


Furnace  Temperature  -  1930°R 


0  5  10  15  20 


Time,  9,  sec 

a.  0.05-in.  Aperture 

Fig.  16  Black-Body  Cavity  Temperature-Time  Histories 
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AEDC-TDR-63-93 


Furnace  Temperature  -  1914°R 

o 

Q  =6.39  Btu/ft  -sec 

source 


b.  0.15-in.  Aperture 
Fig.  16  Concluded 
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AEDC-TDR-63-93 


-  Calculated  Source  Heat  Rate 

Kg  ”  0.035  Btu/f t^-hr-deg/ft 
CpB  =0.22  Btu/lbm-deg 


0  5  10  15  20 

Time,  9,  sec 


a.  0.05-in.  Aperture  (Time  Zero  at  3.75  sec) 

Fig.  17  Black-Body  Cavity  Heat  Rate-Time  Histories 
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Heat  Rate 


AEDC-TDR-63-93 


— —  Calculated  Source  Heat  Rate 
Kg  -  0.04  Btu/f t2-hr-deg/ft 
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•p 

CO 

<0 

2  500 
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400 

300 

200 

100 
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Cpg  =0.23  Btu/lbm-deg 


0  5  10  15  20 

Time,  0,  sec 


b.  0.15-in.  Ap«rture  (Time  Zero  at  2.5  sec) 
Fig.  17  Continued 
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Heat  Rate  x  10  ,  Btu/sec 


AEDC-TDR.63-93 


- Calculated  Sourch  Heat  Rate 

Kg  =  0.04  Btu/f t^-hr-deg/ft 
C  _  =0.23  Btu/lb  -deg 

pB  III 


Time,  0,  sec 


c.  0.15-in.  Aperture  (Time  Zero  at  1.25  sec) 


Fig.  17  Concluded 
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A  E  D  C-T  D  R-63-93 


Qen„_fl0  =6.7  Btu/ft2-sec 
source 


Time,  Q,  sec 

a.  0.05-in.  Aperture 

Fig.  18  Black-Body  Cavity  Temperature-Time  Histories  (External  Convection) 
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Temperature 


1300 


Q  «■  6.4  Btu/ft^-sec 

source 


1100 


900 


700 


500 


Time,  d,  sec 


b.  0.15-in.  Aperture 
Fig.  18  Concluded 


Heat  Rate  x  10  ,  Btu/sec 


AEDC-TDR-63-93 


— Calculated  Source  Heat  Rate 
Kg  ”  0.04  Btu/ft2-hr-deg/ft 
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a.  0.05-in.  Aperture 

Fig.  19  Black-Body  Cavity  Heat  Rate-Time  Histories  (External  Convection) 
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Heat  Rate  x  10' ,  Btu/sec 


AEDC-TDR.63-93 


- —  Calculated  Source  Heat  Rate 

'  Kg  -  0.04  Btu/f t2-hr-deg/ ft 
•  C  B  -  0.23  Btu/lbm-deg 


Time,  0,  sec 

b.  0.15-in.  Aperture 
Fig.  19  Concluded 
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